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Abstract. In this paper we review recent progress in our understanding of the chemical evo- 
lution of protoplanetary disks. Current observational constraints and theoretical modeling on 
' ' ■ the chemical composition of gas and dust in these systems are presented. Strong variations of 
^^ ' temperature, density, high-energy radiation intensities in these disks, both radially and verti- 
A^-^ , cally, result in a peculiar disk chemical structure, where a variety of processes are active. In hot, 
fvj ■ dilute and heavily irradiated atmosphere only the most photostable simple radicals and atoms 
, ' and atomic ions exist, formed by gas-phase processes. Beneath the atmosphere a partly UV- 
r-{ , shielded, warm molecular layer is located, where high-energy radiation drives rich ion-molecule 
h— 5 and radical-radical chemistry, both in the gas phase and on dust surfaces. In a cold, dense, dark 

^vj disk midplane many molecules are frozen out, forming thick icy mantles where surface chemistry 

rvj . is active and where complex polyatomic (organic) species are synthesized. Dynamical processes 
' affect disk chemical composition by enriching it in abundances of complex species produced via 
slow surface processes, which will become detectable with ALMA. 
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2 !• Introduction 



Pro toplanetary dis ks (PPDs) are ubiquitous around young stars (e.g.. lLin fc Papaloizou 



^ • Il980t lLissaueiill987 1. Their chemical composition and physical properties regulate the 



efficiency and timescale of planet formation. Molecules and dust serve as heating and 
T— I cooling agents of the gas, while dust grains also dominate the disk opacities. Molecular 

^ lines observed at infrared and (sub-)millimeter wavelengths are useful probes of physical 

^^ conditions (temperature, density, kinematics, turbulence) in PPDs. 

. -- Over the past decade significant progress has been achieved in our understanding of 

^^ disk chemical composition, both theoretically and observationally. Since molecular hy- 

drogen is not observable when it is cold, we have to rely on other trace species to infer 
disk physics and chemistry (Table [ij. Multi- molecule, multi-transition interferometric 
observations, coupled to line radiative transfer and chemical modeling, allowed to con- 
strain disk sizes, kinematics, dist ribution of tempera ture, surface density, and m olecular 
column densities (see reviews bv lBergin et al.l (|2007l ) and lDutrev et al.l (|2007al )). 



• '^ . Observations of thermal dust emission are employed to constrain dust properties, infer 

1^^ ' disk masses, whereas IR spectroscopy is used to study mineralogical composition of dust 

^ . in PPDs. Recently, with space-borne (Spitzer) and ground-based (Keck, VLT, Subaru) 

■ " " infrared telescopes, molecules have been detected in very inner zones of planet-forming 

systems, at r < 1 — 10 AU. 

The conditions of planets formation in the early Solar system have been revealed by 
a detailed analysis of chemic al and mineral ogical composition of meteoritic samples and 



cometary dust particles (e.g., lBradlevll2005l ) . The recent Stardust and Genesis space mis- 
sions have returned first samples of pristine materials, likely of cometary origin, showing a 
complex structure of high-temperature crystalline silicates embedded in low-temperature 
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Table 1. Molecular species commonly utilized to study disks 



Tracer 


Quantity 


Midplane 


Molecular 
layer 


Atmosphere 


Inner 
zone 


'"CO, "CO 


Temperature 


mm* 


mm 


mm 


IR 


Ha 


— 











IR 


NH3 


— 


cm 


cm 








CS, H2CO 


Density 





mm 





IR 


CCH, HON, CN 


Photochemistry 





mm 





IR 


HCO+ 


Ionization 





mm 








N2H+ 


— 


mm 











C+ 


— 








IR 


IR 


complex organics 


Surface 
processes 


IR** 


IR-cm 





IR,mm 


DCO+, DCN, 


Deuterium 


mm 


mm 








H2D+ 


fractionation 











* - "mm/cm" and "IR" mean radio-interferometric and infrared observations, respectively. 
** - Complex molecules frozen onto dust surfaces could be detected through absorption lines 
in infrared, while the gas-phase counterparts emit at (sub-) millimeter frequencies. 



condensates (JBrownlee et al.ll2004lFlvnn et al.lbooeHBrownlee et al.ll2008l) . The presence 
of crystalline s ilicates in outer regions of protoplanetary disks has also been revealed (e.g., 



van Boekel et al. 2004: Juhasz et al. 2010 ) . An isotopic analysis of refractory condensates 



in unaltered chondritic meteorites shows strong evidence that the inner part of the Solar 
Nebula has been almost completely mixed during the first several Myr of evolution (e.g., 
IBoss 2004: Ciesla 2009) . 

These intriguing fi ndings are partly understood in modern astrochemical models of 
protoplanetar y disks ( Aikawa & Hcrbst 199!?; Markwick et al. 200z; van Zadclhoff e t al.l 



2003: Tscharnuter fc Gail 2007: Agundez et al. 2008: Woods fc Willacy 2009: Visscr et alj 



2010t ISemenov fc Wiebe]|201l[ ). The major result of the chemical 



modeling is that disks have a layered chemical structure due to heavy freeze-out of gas- 
phase molecules in the cold midplane and their photodissociation in the atmosphere. 
Observed column densities of CO, HCO+, N2H+, CN, HON, HNC, CS, etc. are qualita- 
tively agree with the chemical models. 

In this paper we briefly review the major observational findings and modeling predic- 
tions for the chemical composition and evolution of protoplanetary disks. 



2. Observational constrains on composition of gas and ices 

Apart from CO and its isotopologues, and occasionally HCO+, DCO+, CN, HCN 
DCN, CCH, H2 CO, and CS, the molecular content of prot o planetary disks remains largely 
unkriown (e.g.. lOutrev et al.l [l997l: iKastner et al.l Il997l: lAikawa et al.l l2003t iThi et al 



2004 IPietu et al.ll2007tlOi et al.ll2008t iHenning et al.ll2010t ). Molecular line data are lim 
ited in sensitivity and resoluti on. Thus the spatial distribution of molecular abundance s 
is still poorly determined (e.g.. |Pietu et al.ll2005l : lOutrev et"aLll2007bl : [Panic et al.ll2009l ). 
Observational facilities such as the Plateau de Bure interferometer (France) and Sub- 
millimeter Array (USA) have permitted measurements of several brightest nearby disks 
(DM Tau, LkCa 15, AB Aur, and TW Hya). Soon, with the Atacama Large Millimeter 
Array entering it's full power in 2013, we will make substantial progress in our knowledge 
of molecular disk structure. 

Typically, studies of disk physics begin with observations of bright CO lines. These lines 
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are thermally excited at densities ~ 10—10 cm •^. The CO lines are optic ally thick and 



their intensities measure kinetic temperature in the upper disk layer (e.g.. lDutrev et alJ 
19971). The lines of less abundant ^'^CO and C^^O are typically optically thin or partially 
optically thick and are sensitive to both temperature and corresponding column densities 
throughout the entire disk. Strong CO lines are most suitable for accurate determination 
of disk kinematics as well as orientation and geometry. Their measured widths indicate 
that turbulence in d i sks is subsonic, with t ypical velocities of about 0.05 — 0.2 kms"""^ 



( Dutrev et al.ll2007al: ISemenov et al.ll2010al: iHughes et al.ll2011[ ) 



It has been found that disks appear progressively larger from observations of dust 
continuum, onwards to C^^O, ^"^CO, and ^^CO, respectively, with typical values of 
300 — 1000 AU. This is a manifestation of selective isotopic photodissociation. Most 
disks exhibit a vertical temperature gradient, ranging from '^-^ 10 K at the midplane 
to ^ 50 — 100 K i n the a tmosph ere region, as determined by physi cal models (e.g., 
Dartoiset all 120031: lOi et al . 2006; Pietu et all 120071: llsella et al.ll2010[ ). though several 
disks with lar ge inner cavities do not show evidence for such a gradient (e.g., CM Aur 
and LkCa 15: lDutrev et al.ll2008l: iHughes et al.ll2009l ). 

By comparing intensities of 6-5 to 2-1 CO transitions it was found bv lOi et al.l ( 20061 ) 
that the TW Hya disk has a surface region that is superheated , so that an additiona l 



heating mechanism is required, possibly stellar X-ray radiation ([Glassgold et al.l 119971 ) 



A significant reservoir of very cold CO, HCO+, CN and HCN gases has been found in the 
disk of DM Tau at temperatures < 6-17 K, which cannot be explained by conventional 
chem ical models w ithout invoking a non-thermal dcsorption or transport mechanism 
(e.g.. 



Semenov et al. 2006: Aikawa 2007 : Hcrsant et al. 2009) . 



The second most readily observed molecular species in disks is HCO"'". The low- lying 
transitions of this ion are thermalized at densities of about 10^ cm~^. This is one of the 
most abundant charged molecule in disks, the other being C^ (not observable at millime- 
ter wavelengths) and Hg" (lacking dipole moment). The ioni zation degree measured from 
HCO+ is ^ 10~^° — 10~^ inside warm molecular disk layer (|0i et al.lbOOSt iDutrev et al. 
2007bl) . 

Less strong lines of C2H, CN, and HCN are sensitive to the intensity and shape of 



the inc ident UV spectrum and are excellent tracers of photochemistry (e.g.. lBergin et al 



2003al ). Recently detected DC O"*" and DCN h ave abundances that are about 1 — 10% of 
the HCO"*" and HCN densities ( Oi et al.ll2008l ). It remains to be verified whether such a 
large degree of deuteration is a heritage of cloud chemistry or produced in situ. 

A key observational result is that molecular abundances are deplet ed by factors 5- 
100 compared to the values in the Taurus Molecular Cloud (JDutrev et al. 20073,). Since 
disks have higher densities up to 10^ — 10^" cni~^ and enshrouded in stronger ionizing 
radiation fields, this depletion can be attributed to a combined effect of photodissociation 
and freeze-out. 

The results from the Infrared Space Observatory and Spitzer telescope have proven the 
existence of a significant amount of frozen mat erial and various types of silicates and poly- 
cyclic aromatic hydrocarbon s (PAH) in disks ( van den Ancker et al.ll2000l : Ivan Dishoeck 



2004t iBouwman et al.ll2008l) . The PAH features at ~ 3-12/i m probe the inciden t ra- 



diation field and density distribution of the upper disk (e.g., iHabart et al.l 120041 ). Re 
cently, with space-borne (Spitzer) and ground-based (Keck, VLT, Subaru) infrared tele- 
scopes, molecules have been detected in very inner zones of planet-forming systems, at 
r < 1 — 10 AU. Rotational-vibrational emission lines from CO, CO2, C2H2, HCN, OH, 
H2O imply a rich chemistry drive n by endothermic reactions or reactions with activa- 
tion barriers and photoprocesses ( Lahuis et al.l l2006t ICarr fc Naiitall2008t ISalvk et al 
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Figure 1. Physical and ciiemical structure of a protoplanetary disk. Timescales are rough 
estimates and given for a radius of about 100 AU. 



2008t|Pontoppidan et alll2008HPascucci et al.|[2009tlvan der Plas et al.|[20oHSalvk et al 



201 ih . Through ISO and Spitzer infrared spectroscopy abundant ices in cold disk regions 
consisting of water ice and substantial amounts (~ 1 — 30%) of vo latile materials like 
CO CO2, NHcj, CH4. H9CO, and HCOOH ha ve been detected (e.g.. |Pontoppidan et aT 
2005HTerada et al.ll2007t IZasowski et al.ll2009 i). 



In Table [T] the various molecules used to study protoplanetary disks are summarized. 



3. Global evolutionary picture of disk chemical structure 

The current scheme of the disk chemical structure is presented in Fig [1] The disk 
can be divided into 4 distinct chemical regions, primarily determined by temperature. 
The warm "inner zone" corresponds to radii < 20 AU that are accessible with the IR 
instruments, while the other 3 regions dubbed "midplane", "molecular layer", and "at- 
mosphere" represent the outer disk regions (r > 20 — 50 AU) that are observed with 
radio-interferometers . 

Dense midplane is opaque to the UV and X-ray radiation and remains cold and essen- 
tially neutral (10 — 20 K). Chemical complexity in this region is initially reached by fast 
ion-molecular reactions in the gas, followed by slow accretion onto grains and surface 
reactions (mostly hydrogenation) . Once formed, a molecule is seldom re-emitted into the 
gas phase. Typical chemical timescales for this region are determined by freeze-out and 
surface reactions, and are > 10^ — 10^ years in the outer disk. In the very inner part, the 
midplane is hot due to viscous accretion heating (T > 50 K) and gas-grain interactions 
are not important, so the chemical timescale is due to gas-phase neutral-neutral reactions 
(~100 years). The planet-forming zone is in general in a chemical steady-state, which 
allows the application of simple condensation/evaporation thermodynamical chemical 
models. 

Adjacent to the midplane lies a less dense, warmer layer (T ~ 30 — 70 K) that is partly 
shielded from stellar and interstellar UV/X- ray radiation. Young stars emit intense non- 
thermal UV and thermal and variable X-ray radiation field s that strongly mediate disk 
chemistry at intermediate heights fe.g.. iBergin et al.l 20071). UV e xcess in T Tauri stars 
is related to their chromospheric activity fe.g.. lBouvier et al.ll2007l ). Soft X-ray radiation 
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Table 2. Chemical reactions active in disks 

Process Formula Midplane Molecular Atmosphere Inner 

layer zone 

Bond formation 

Radiative association A + B — >■ AB + hf X X XX 

Surface formation A + B||gr -)■ AB + gr X X 

Three-body A + B + M^AB + M X 



Bond destruction 












Photodissociation 


AB + hi^ ^ A + B 





X 


X 


X 


Dissociation by CRP 


AB + CRP ^ A + B 


X 


X 








Dissociation by X-rays 


— 





X 


X 


X 


Dissociative 


AB+ -1- e" ^ A + B 


X 


X 


X 


X 


recombination 












Bond restructuring 












Neutral-neutral 


A + BC ^ AB + C 


X 


X 





X 


Ion-molecule 


A+ + BC ^ AB+ + C 


X 


X 


X 


X 


Charge transfer 


A+ + BC -> A -h BC+ 


X 


X 


X 


X 


Unchanged bond 












Photoionization 


AB + hiy -^ AB+ -1- e" 





X 


X 


X 


Ionization by CRP 


AB + CRP ^ AB+ + e- 


X 


X 








Ionization by X-rays 


— 





X 


X 


X 



has likely the same origi n, while ha rder X-rays are produced in magnetic reconnection 



loops or due to jets (e.g.. lGiidel fc N aze 2009) 



These energetic photons dissociate and ionize gas and photodesorb surface species, 
thus enriching the gas composition and initiating a rich chemistry. Abundances of most 
molecules attain high concentrations in this zone, and numerous molecular lines are 
excited and observable emission produced. Chemistry does not reach a steady state in 
this region and a typical timescale, as determined by surface chemistry and gas-grain 
interactions, exceeds 10^ years. 

Above the intermediate layer a hot, dilute, and heavily irradiated disk atmosphere 
exists (T > 100 K). This is a molecularly-poor region, where only simple light hydrocar- 
bons, their ions, and other radicals such as CCH and CN are able to survive. Chemical 
timescales are short (^100 years) and defined by photochemical processes and dissociate 
recombination. 



4. Gas-phase chemistry 

Models of chemical evolution employ a set of chemical reactions of relevance to con- 
ditions in disks. Modern astrocherri i cal databases include up to 600 species i r ivolved in 
4-700 reactions fem ith et all 120041: lie Petit et all 120061 IWoodall et alll2007t IWakelam 



2009l) . Only 10 — 20% of reaction rates have been studied in the laboratory or cal- 



culated theoretically and thus models are prone to inherent uncertainties (see, e.g., 
Wakelam et al.ll2010l ). All reactions can be divi ded into 4 distinct g roups that are domi- 



nant in different disk regions, see Table [2] (after Ivan Dishoeck &: Black 1988) . 

Apart from the very dense inner zone, all reactions in disks are two-body processes. 
Three-bod y reactions become competitive only at < 10 AU, where n > 10^° cm~'^ 
( Aikawa et al. 19991 ). The main processes leading to formation of molecular bonds are 



slow radiative association and surface reactions. Upon collision, a coUisional complex 
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in an excited state ma; 
of a photon (e.g., 



form, which is stab i lized with a low probabili ty by emission 

19721 : iHerbst fc KlempereiHToTl . For example, 

formation of light hydrocarbons starts with radiative association of C"'" and H2, leading 



e may rorn i , wnicn is 
BateslllQSlI: [Wilhamd 



to excited CH+ (,Hcrbstlll985[ ). 



Ionized by FUV, X-ray photons or CRPs, ions and molecule drive rapid ion-neutral 
chemistry, which constitutes the largest fraction of astrochemical models. These reactions 
are exothermic, with high rate co efficients ~ 10~^ cm~^ s~^, which often increase toward 



low temperatures {(3 < 0) fe.g.. IPalgarno fc Blackl 119761 ). fon-neutral reactions result 
in bond restructuring of the reactants. Some of the most important reactions of this 
category are protonation reactions. 

Molecular ions are destroyed by dissociative recombination with electrons and neg- 
atively charged grains. These processes are es pecially fast at low temperatures, with 



typical rates of about 10 "cm ^s ^atlOK (jWoodall et all 120071 ). For nearly all ob 



served species, dissociative recombination is an important formation pathway (e.g., water 
and hydrocarbons). Often, at later evolutionary times, > 10^ years, dissociative recom- 
bination is balanced by protonation reactions, e.g. CO -I- H^^ — > HCO+ -|- H2 followed 
by HCO+ + e^ — > CO -I- H. Products and branching ratios for pol yatomic ions are not 
easily obtainable fe.g.. lBates fc Herbstlll988l:ISpanel fc Smithlll994l) . 

A number of neutral-neutral reactions involvin g open-shell radicals can also be active 
in both cold outer and warm inner disk regions ( van Dishoecklll998l ). The typical rate 



coefficient for these reactions is ^ 10~ -10~ cm~ s~ \ i.e., only about an order of 
magn itude lower than for the ion-molecule processes (e.g. JClarv et al.lll994ilSmith et al. 
20041 ). One of the most interesting reactions of this ty pe is formation of H CO+ upon 



collision between O and CH (ao = 2.0 • 10"" cm""^ s~^: IWoodall et al.ll2007l ) 



5. Photochemistry 

The photochemical reactions are overviewed in this book in the paper by van Dishoeck 
& Visser and only briefly discussed here. 

Young T Tau stars {Tcs — 4 000 K) emit inte nse non-thermal UV radiation that has 
a spectra different from the interstellar UV field (JBergin et al.ll2003al ) , with a prominent 



Lymaua line, while hot (Teff > 10 000 K) Herbig Ae/Be stars produce a lot of thermal 
UV emission. The overall intensity of the stellar UV radiation at 100 AU from the star 
can be as high as 500 and 10f_Jor aT Tau and a Herbig Ae star, respectively, in units of 
the interstefiar UV held (|Habing|ll968l ). 

Many molecules, such as CO, H2, and CN, are dissociated by radiation at short 
wavelengths (A < 1 100 A) , while photodissociation of other species, such as HC N, oc- 
curs at longer wavelengths ([van Dishoeck fc Blacklll988t Ivan Dishoeck et al.ll2006l) . Some 
molecules (e.g., H2 and CO) dissociate via absorption of UV photons at discrete lines 
to the excited (Rydberg) states, whereas other molecules are disso ciated either by the 
continuum (e.g., CH4) o r by the continuum and fines (e.g., C2) (jvan Dishoeckl 1 1988 : 
van Dishoeck et al.ll2006l ). The high ratio of CN to HCN abundances, observed in disks 



(jDutrev et al.lll997l ). can be explained by mor e intense photodisso ciation of HCN if part 
of the UV flux comes as Lyo, photons fl216A: lBergin et al.ll2003al ). 

Since the dissociation of the abundant H2 and CO molecules result from photoab- 
sorption at discrete wavel engths, isotopically selective photodissociation based on self- 
shielding is possible (e.g.. lVisser et al.ll2009al ). Two conditions for that are required: 1) 



dissociation via line absorption for each isotopically substituted molecule, and 2) differ- 
ential photolysis that depends upon the isotopic abundances. Self-shielding occurs when 
the spectral lines leading to dissociation of the major isotopic species optically saturate, 
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while the other residual lines relevant for dissociation of the minor isotopes remain trans 
parent. Under conditions of proto planetary disks, isotopic-sele c tive chemistry may occur 
for abundant CO molecules (e.g., 
1983t lyan Dishoeck fc Blacklll988 : 



Dalgarno fc Stephens 1970; Thie mens fc Heidenreich 
Lee et al.lll996; .Visser et al...2od9a[K 



6. Gas-grain interactions 

In cold disk regions grains serve as a passive sink for heavy molecules and electrons, 
while providing free electrons to dissociate positi ve ions. In the darke st midplane zone 
dust grains become the dominant charged species ( Semenov et al.ll2004[ ). At low temper- 
atures of < 20 — 50 K, many molecules adhere to a grain with a nearly 100% sticking 
probability since their kinetic energy is much smaller than the binding energy (see, e.g., 
d'Hendecourt et al.lll985HBuch fc Zhane]|l991 ). This sticking probability also depends on 
surface properties, such as porosity and distribution of surface chemi- and physisorption 
sites. Chemisorption requires formation of a chemical bond between a surface sp ecies and 
a grain, and such a species does not evaporate easily (e.g.. lCazaux et al.ll2005l ). 

The most effective desorption processes for disk chemistry are thermal evaporation, 
cosmic ray induced desorption, and photodesorption. Thermal evaporation occurs if a 
molecule has energy that exceeds it's binding energy. Typical binding energies of ph- 



ysiso rbed species are about 1 000 K for light molecules like CO and N2 (jBisschop et al 



20061 ) ■ and much larger for heavier cyanopolyynes and carbon chains. Chemisorbed 
species do not desorb until very high temperatures of 100 — 1000 K are reached. 

In disk midplane CRPs provide energy that allows dust mantles (partly) thermally 
evaporate. A relativistic iron nucleus may eventually collide wi th a grain and impu l sively 
heat it, releasing a portion of the volatile co mponent (e.g.. IWatson fc Salpeteii Il972l : 
Leger et al.lll985t iHartauist fc Williamslll990l) . In less opaque disk regions, penetrating 
UV photons lead to photoevaporation of surface species. The probability of evaporation 
per one UV photon has been measured in the laboratory for some simple molecules , 
such as CO, H2O, CH4, and NH3, and is about 10-6-10^2 (e.g.. lOberg et all 12007 . 
2009blla). The dilute UV radiation field produced by cosmic rays in the disk midplane 
(jPrasad fc Tarafdaiill983f ) can also be important for high values of the photodesorption 
rates. Another me chanism is the X-ra y induced desorption which is efficient only for tiny 
grains, < 0.05/im (JNaiita et al.ll200l[ ). 



7. Surface formation of complex species 

The surfaces of dust grains serve as a catalyst for reactions that do not proceed ef- 
ficiently in the gas phase. The most notable example is f ormation of molecular hydro 



gen, which occurs almos t entirely on dust surfac es (e.g.. iHollenbach fc Salpeteii 11971 : 
Watson fc Salpeteiill972l ). iTielens fc Hagenl ()l982h have studied the chemical evolution 
of surface species on large grains and considered the effect of various migration rates for 
light atoms and heavy molecules. A typical O.lfim amorphous silicate grain accommo- 
dates ^ 10^ surface sites available for accretion. An atom or light radical, if it is not 
chemisorbed, may migrate over the surface from site to site by thermal hopping, when 
its energy exceeds the barriers for particle motions and react with other species. 

There is an increasing body of evidence from laboratory measurements that com- 
plex molecules like methanol cannot be produced in the gas phase via radiative asso- 
ciation forming a lar ge protonated precursor followed by its dissociative recombination 
( Geppert et al.ll2005[) . Thus, surface reactions due to thermal hopping remain the only vi- 
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able formation pathway for production of complex (organic) molecules in protoplanetary 
disks, found also in meteoritic samples in Solar system. 



8. Importance of dynamical processes for disk chemistry 

While most of the chemical studies are still based on laminar disk models, evidences 
for mixing call for a more sophisticated treatment. Models of t he early Solar nebula 
with radia l transport by advective flows have been developed (e.R.,lMorfill k. Volklll984l: 



Gail 120011 [JnoS IWehrstedt fc GaillEonS iBossllJOoS iKeller fc Gailll2nn4[ ). Illgner et al 
( 20041 ) for the first time modeled the influence of turbulent diffusion in the vertical 



direction and advection flows in the radial direction on the chemical composition of the 
inner disk region. They found that dynamical processes s ignificantly affect the chemical 
evolution of sulfur-bearing species. IWillacv et al.l (20061) have shown that ID vertical 
mixing modifies chemical com position of the o u ter di sk region and t hat the m i xing r esults 
better agree to observations. ISemenov et al.l ( 20061 Paper I) and lAikawal ( 2007 ) have 
found that turbulent transport allows expla ining the presence o f a larg e amount of cold 
(< 15 K) CO gas in the disk of DM Tau. iTscharnuter fc Gaill (|2007t ) have used a 2D 
disk chemo-hydrodynamical model and showed that in the disk midplane matter moves 
outward, carrying out the angular momentum, while the accretion flows toward the star 
are located at elevated altitudes. Consequently, gas-phase species produced by warm 
chemistry in the inner nebula can be steadily transported into the cold outer regio n 
and freeze out. A radial advection model has also been utilized bv lNomura et al.l ( 2009I ). 
who have demonstrated that inward radial transport enhances a bundances of organic 
molecules (produced mainly on dust surfaces in cold outer regions). iHersant et al.l ( 20091 ) 
have studied various mechanisms to retain gas-phase CO in very cold disk regions. They 
concluded that efficient photodesorption in moderately obscured disk regions (Ay < 5'") 
greatly enhances g as-phase CO concentra tions, while the role of vertical mixing is less 
important. Finally, iHeinzeller et al.l (|201l[ ) have investigated the disk chemical evolution 
with radial advection, vertical mixing, and vertical wind transport processes. They have 
found that the disk wind has a negligible effect on disk chemistry, whereas the radial 
accretion alters the molecular abundances in the cold midplane, and the vertical turbulent 
mixing affects the chemistry in the warm molecular layer. 

We study t he influence of turbulence mixing on the chemical evolution of protoplan- 
etary disks in ISemenov fc Wiebd (2011). Mixing is important in disks since a chemical 
steady-state is not reached due to long timescales associated with surface chemical pro- 
cesses and slow evaporation of heavy molecules. Our analysis was based on the a-model 
of a ~ 5 Myr DM Tau disk cou pled to the large-scale gas-grain chemical code "AL- 
CHEMIC" ([Semenov et al.ll2010bl ). To account for production of complex molecules, our 
chemical network was supplied with a set of surface reactions and photoprocessing of 
ices. 

The adopted fla r ing d isk structure is based on a 1-l-lD steady-state a-model similar 
to lD'Alessio et al.l (j 199 9 1 model. T he non-therma l FUV radiation field from DM Tau is 
represented by the scaled ISRF of Draine (llQTST). w ith the un-attenuated intensity at 
100 AU of x*(100) = 410 (e.g., iBergin et al.ll2003b[) . For the X-ray luminosity of the 



star we adopt a value of 10 ergs , which is constrained by recent measurements with 
Chandra and XMM in the range of 0.3-10 keV. 

The turbulence in disks is likely driven by the magnet orotational instability ( MRI), 
which is operative even in a weakly ionized medium fe.g.. lBalbus fc Hawlevlll99ll ). This 
turbulence causes anomalous viscosity that enables e fficient redistribution of the angular 
momentum. We have followed the parametrization of IShakura fc SunvaevI (|l973l ). where 
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Table 3. Detectable tracers of turbulent mixine. 



Steadfast 


Hypersensitive 


CO 


Heavy hydrocarbons (e.g.jCsHe) 


H2O ice 


C2S 




C3S 




CO2 




O2 




SO 




SO2 




OCN 




Complex organics (e.g., HCOOH) 



turbulent viscosity 1^ is related to local disk properties such as the vertical spatial scale 
H{r), the sound speed Cs{r,z), and the dimensionless parameter a: 

v{r,z) = acs{r,z) H(r). (8-1) 

From observational constraints a is ^^ 0.001—0.1 ( Andrews fc Williamsll2007tlGuilloteau et al 



20111) . so we adopt the constant value of 0.01. Consequently, the diffusion coefficient is 



calculated as 

Aurb(r,^) = z^(r,z)/5c, (8.2) 

where Sc is the Schmidt number describing efficien cy of turbulent diffusivity (see e.g. 
Shakura k. Sunvaevlll973l : ISchrapler fc Henninej|2004l) . We assume that gas-phase species 



and dust grains are well mixed, and transported with the same diffusion coefficient. 

We show that the higher the ratio of the characteristic chemical timescale to the turbu- 
lent transport timescale for a given molecule, the higher the probability that its column 
density will be affected by dynamical processes. Thus, turbulent transport influences 
abundances of many gas-phase species and especially ices. Vertical mixing is more im- 
portant as it affects the evolution of gas-phase and surface species of any kind, whereas 
the effect of radial mixing is pronounced mostly for the evolution of ices. The radial 
temperature gradient is weaker, and thus is only relevant for the evolution of polyatomic 
ices formed via surface reactions of heavy radicals, whereas steep vertical gradients of 
temperature and high-energy radiation intensity cause much sharper transition from the 
ice-dominated chemistry in the disk midplane to the rich gas-phase chemistry in the 
molecular layer. 

The simple molecules that are unresponsive to dynamical transport include C2H, C+, 
CH4, CN, CO, HCN, HNC, H2CO, OH, as well as water and ammonia ices. The sensi- 
tive species to dynamics are carbon chains and other heavy species, in particular sulfur- 
bearing and complex organic molecules frozen onto the dust grains. Mixing steadily 
transports ice-coated grains in warmer regions, allowing more efficient surface process- 
ing due to enhanced hopping rates of heavy radicals. In warm intermediate layer these 
organically-rich ices evaporate, and in the inner disk they can also be photodissociated 
by CRP/X-ray-induced UV photons. The importance of mixing is higher in an inner, 
planet-forming disk zone, where thermal, density, and high-energy radiation gradients 
are stronger than in the outer region. Still, mixing does not completely erase the layered 
chemical structure of protoplanetary disks. 

Several promising detectable tracers of dynamical processes in protoplanetary disks 
are the column density ratios of the CO2, O2, SO, SO2, C2S, C3S to that of CO and the 
water ice (see Table[3]). The detection of complex species (e.g., dimethyl ether, formic acid. 
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methyl formate, etc.) in protoplanetary disks with ALMA and JWST will be a strong 
indication that chemical evolution of these objects is influenced by transport processes. 
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Discussion 

Paul M. Woods: There were some dynamical models recently that show that the accre- 
tion flow proceeds inwards along the disc surface and outwards along the disc midplane. 
Are those models credible, and how would that affect the disc chemistry? 

D. Semenov: Turbulence is a 3D phenomenon and its magnitude varies greatly across 
a disk. Global 3D MHD models show that viscous transport occurs in all directions in 
any given disk location, albeit with different efficiencies, and that there is no a specific 
direction for it. According to observations, advective transport seems to be an effective 
mechanism that regulates overall lifetime of disks, ^1 — 10 Myr, in which most of disk 
matter accretes onto the central star and photoc vapor ates. Consequently, disk structure 
changes with time, as well as the FUV and X-ray penetration, and grain properties. 
This strongly affects chemical evolution of protoplanetary disks, and favors the use of 
evolutionary disk models (see R. Vissier's presentation and proceeding paper. 

James R. Lyons: What turbulent viscosity parameter a-values were used for slow vs. 
fast mixing? 

D. Semenov: We considered Fickian diffusion, in which the diffusion coefficient -Dturb(?', z) 
is attributed to the a-parameterized viscosity, a = 0.01: Z?turb(?'jz) = h'{r,z)/Sc = 
aCs{r, z) H{r)/Sc, where the characteristic spatial scale is H{r) and the sound speed is 
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Cs{r,z). The Sc the Schmidt number describing efficiency of turbulent diffusivity. The 
"fast" mixing model has Sc^l and the "slow" model has "Sc=:100". 

James R. Lyons: Have you included isotopes in vertical mixing calculations? 

D. Semenov: No. 

James R. Lyons: What did you mean by "isotopic homogeneity of the Solar nebula"? 
It is not homogeneous in oxygen isotopes. 

D. Semenov: I've meant that the inner, 1-20 AU region of the Solar nebula has been 
isotopically homogenized at bulk level, though oxygen anomalies at percentage level still 
persist. 



